Interactions Between Rydberg-Dressed Atoms 



J. E. Johnson, S. L. Rolston 

Joint Quantum Institute, Dept. of Physics, University of Maryland and National 
Institute of Standards and Technology, College Park, MD 20742-4111, U.S.A. 

(Dated: August 24, 2010) 

We examine interactions between atoms continuously and coherently driven between the ground 
state and a Rydberg state, producing "Rydberg-dressed atoms." Because of the large dipolar cou- 
pling between two Rydberg atoms, a small admixture of Rydberg character into a ground state can 
produce an atom with a dipole moment of a few Debye, the appropriate size to observe interesting 
dipolar physics effects in cold atom systems. We have calculated the interaction energies for atoms 
that interact via the dipole-dipole interaction and find that due to blockade effects, the independent 
two-atom interaction terms are limited in size, and can be ii-independent up until the dipolar energy 
is equal to the detuning. This produces independent interactions different from the expected 1/R 3 
dipolar form, which have no direct analogy in condensed matter physics, and could lead to new 
quantum phases in trapped Rydberg systems. 



I. INTRODUCTION 

There have been numerous theoretical predictions of 
novel phases of matter in ultracold atomic systems with 
long-range interactions, including dipolar crystals [UHJ, 
supersolids, [IHS], striped and checkerboard phases [5]- 
|7], and others [5]. These calculations have generally as- 
sumed the existence of ground-state polar molecules with 
dipole moments in the range of 2 — 5eao- Although there 
have been dramatic accomplishments recently in the pro- 
duction of ground-state polar molecules [Bj, they have 
yet to be used for such dipolar many-body physics, and 
are challenging to produce. As has been pointed out 
by [I], Rydberg atoms might also be able to fill this role. 
The maximum dipole-dipole interaction between two Ry- 
dberg atoms with principle quantum number n is of or- 
der n 4 eao, orders of magnitude larger than needed for 
predicted dipolar effects, and in fact so large that the in- 
terparticle forces would overwhelm any optical trapping 
forces from an optical lattice, for example. In addition, 
the typical lifetime of a Rydberg state with n = 50 is 
about 100 /is too short to allow a many-body sys- 
tem to reach equilibrium. Because the full Rydberg- 
Rydberg interactions are so strong, we can use a state 
with only a fraction of those interactions, something we 
can achieve by creating a wavefunction that is mostly 
ground state with a small, adjustable Rydberg compo- 
nent, using a coherent coupling, dressing the atom |12j . 
This would be accomplished with continuous laser irra- 
diation of ground state atoms, coherently coupling the 
ground state to the Rydberg state via a one- or two- 
photon transition. We can imagine creating atoms that 
have only 1% Rydberg character, which is still sufficient 
to create interesting dipolar physics, and increases the 
lifetimes to ~ 10 ms, where it may be possible for the 
system to come to equilibrium. The admixture fraction 
is controllable by adjusting the coupling laser detuning 
and intensity, and thus would also give a tunable dipolar 
coupling between atoms, which could be a useful fea- 
ture in exploring the effects of long-range interactions. 
Additionally, the dipole-dipole interaction could be de- 



pendent upon an externally-applied static electric field, 
either through tuning near a Forster resonance |10j . or 
by inducing a dipole moment via the Stark Effect [llj . 

In what follows, we will assume an idealized Rydberg 
coupling laser, described solely by a Rabi frequency £1 
and a detuning S with respect to the one-atom Rydberg 
transition. In practice this would be created through a 
two-photon excitation with a large intermediate state de- 
tuning to assure coherent coupling. In |13j they demon- 
strate coherent coupling between a ground and Rydberg 
level in Rb with a two-photon Rabi frequency of ~ 100 
kHz, and an intermediate state detuning of ~ 500 MHz. 
Our goal is to create Rydberg-dressed atoms with a wave- 
function 

\fp)=a\g)+0\r) (1) 

where \g) is the ground state and \r) is the Rydberg state. 
Such a state would have a spontaneous decay rate of 

T~Kff|^>| 2 ~/9 2 7r (2) 

where j r is the Rydberg decay rate, and d is the dipole 
operator for spontaneous emission from the Rydberg 
state [IT]. It is tempting to then simply calculate the 
dipole-dipole interaction between the two dressed states 
as 

tint = (ip\U dd \i>) = p 2 (r\U M \r) = ft r (3) 

where U dd is the usual dipole-dipole operator between 
Rydberg states, and e r is the full interaction energy be- 
tween two Rydbergs, which can be of order 10 GHz for 
R = 1 fim and n ~ 50. As we will see below, this ex- 
pression is in general invalid, because the atom-atom in- 
teractions will cause a blockade effect [14Tll6| . such that 
the two-atom wavefunction contains much less than j3 2 of 
the \r)\r) state. The correct procedure is to calculate the 
dressed states for two atoms simultaneously. The block- 
ade effect will arise naturally out of the dressed eigen- 
states of this two-atom system. 
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FIG. 1: The energy levels of the two-atom state for 

Q/2TT = 10 MHz, (5/27T = 100 MHz, and a.) 
c 3 /2n = -1000 MHzx^rn 3 , and b.) c 3 /2n = 1000 
MHzxfim 3 . 



II. TWO INTERACTING DRESSED ATOMS 

The simplest case of interacting dressed atoms to con- 
sider is the two-atom ease. We could write a 4 x 4 Harnil- 
tonian matrix in the basis \gg), \gr), \rg), |rr), but if in- 
stead we use a basis with l/^/2{\gr) ± \rg)), the anti- 
symmetric state is uncoupled, and can be ignored. The 
resulting 3x3 Hamiltonian in the basis \gg), 1/V2(\gr) + 
\rg)), \rr) can be written as 



H 



( o 







Q 

n_ r n_ 
^ 2S + U dd 



(4) 



If Udd = 0, we recover the dressed states and eigenen- 
ergies for two independent dressed atoms, as expected. 
After diagonalizing the matrix, we find the energy lev- 
els as shown in Fig. [l] where we have used a sim- 
ple C3/R 3 for the dipole-dipole interaction energy, with 
c 3 /2ir = ±1 GHzx/im 3 , typical for a dipole-dipole inter- 
action for n ~ 40 — 50, assuming an applied static electric 
field, for example. We subtract off the constant light shift 
due to the coupling laser so that we plot the interaction 
energy solely due to the two-atom effects (we also ignore 
the angular dependence of the dipole interaction, effec- 
tively assuming a fixed interatomic axis). Figure [I] shows 
two cases, for positive and negative c 3 . In the case of 
negative c 3 there is an avoided curve crossing when the 
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FIG. 2: A detail of the interaction energies for the cases 
shown in Fig. [T] Part b.) shows the convergence with 
an interaction ~ 1/r 3 at longer distances. 



laser is two-photon resonant with the dipole-shifted \rr) 
state, i.e. 25 = —c 3 /R 3 . We consider both the positive 
and negative c 3 cases with the detuning chosen so that 
the interaction energy at short range is positive. Be- 
cause the trapping potential confines the atomic gas, an 
attractive interaction would lead to collapse of the cloud 
rather than emergent order |17j . Figure [2] shows a detail 
of the interaction energies. Note that inside this avoided 
crossing the eigenenergy of the state that connects to the 
ground state (the state of interest in the context of a Ry- 
dberg dressed atom) is almost independent of R, while 
outside the crossing it falls off rapidly with R. For pos- 
itive c 3 (Figs. [I]d, [2}d) there is no avoided crossing, yet 
the eigenenergy also becomes independent of R at dis- 
tances less than about 0.5 /im, which is comparable to 
small optical lattice spacings. This i?-independence is 
a consequence of the blockade phenomenon. The large 
U d d term in the Hamiltonian makes the coupling from 
the single-atom excited state to the doubly excited state 
off- resonant, significantly reducing the \rr) component of 
the two-atom wavefunction (to much less than /3 2 ). 



III. TWO-ATOM EIGENSTATES 

We can find the eigenstates of this matrix, as a function 
of R. Examination of these eigenvectors shows that for 
the eigenvector associated with the lowest energy eigen- 
value, the system is predominantly in the ground state, 
with small amounts of the two excited states mixed in. 
Using these state populations, we can calculate the per- 
centage of Rydberg in the admixed state from the ex- 
pression 



2|<rr#>| 2 



2 1 (5^1 V') 



(rg\ip)\ 



(5) 



which shows that, for the values used previously, the state 
is 0.5% in the Rydberg state, which yields a lifetime of 
about 25 ms. Inside the blockade radius, the system can- 
not be in the state where both atoms are excited to the 
Rydberg level, due to blockade, so the Rydberg charac- 
teristic inside the blockade radius comes from the eigen- 
state where only one of the atoms is excited to the Ryd- 
berg level. 

In the limit when Udd 3> ^ a t small R we can expand 
the ground-state energy of this matrix near R = 0. The 
ground state energy inside the blockade radius becomes 
Egg/k w (l/2)(5-V2ft 2 +8 2 ). Taking the difference be- 
tween this energy and the non-interacting ground state 
energy gives an energy due to interactions of Ei nt /h ~ 
(l/8)(0 4 /<5 3 ), for i > fl. For experiments with cold 
atoms, interaction energies need to be of the order of 
kHz, of the same order as ground state atom-atom in- 
teractions for atoms trapped in a single lattice site (such 
as is relevant for the Bose-Hubbard Hamiltonian). From 
this expression, we can see that it is possible to achieve 
such magnitude of interaction energy with Rabi frequen- 
cies of order 10 MHz and detunings of order 100 MHz, 
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FIG. 3: The energy of the dressed state, for two 
different interactions: the dipole-dipole that varies as 
1/i? 3 (solid), and the Van der Waals that varies as 1/i? 6 
(dashed), for S/2n = 100 MHz, c 3 /2tt = 1000 
MHzx/im 3 , and c 6 /2tt = 500 MHz x fim 6 . 



while at the same time keeping the Rydberg fraction to 
1% or less (necessary for sufficiently small spontaneous 
emission rates). The fi 4 dependence of the interaction 
energy puts a premium on large Rabi frequencies, which 
will be the primary experimental challenge, especially if 
it needs to be implemented over a large spatial volume. 

A second limiting case we can consider is when the 
atoms are separated by a distance much greater than 
the blockade radius, or 26 Udd- We find that the 
ground-state energy eigenvalue dependent on R scales as 
C/dd(l/16)^ 4 /^ 4 . This is as expected: if we can ignore 
blockade, the admixture ratio of the excited state for a 
single atom goes as j3 ~ (1/2)0/(5, so the doubly excited 
state mixture would go as /3 2 , and the matrix element 
(rr\Vdd\rr) scales with /3 4 , as found. Note that once again 
the interaction energy scales as 51 4 . 

The above discussions are independent of the actual 
nature of the atom-atom interaction. Although presented 
in terms of a pure dipole-dipole interaction, ~ C3/R 3 , 
the same blockade mechanism applies for any interaction 
larger than the detuning. Figure [3] plots the energy of 
the dressed state connected to the ground state for two 
interactions : pure dipole-dipole (c 3 /R 3 ) and pure Van 
der Waals (cq/R 6 ). A potential such as that shown in 
[TU] for two interacting 50s Rb atoms that is dipolar at 
short range and van der Waals at long range would fall 
in between the curves shown in Fig. [3] In all cases the 
eigenenergies evolve similarly with i?, with the only dif- 
ference the rate of change with R. 



IV. POSSIBLE EXPERIMENTAL 
IMPLICATIONS 

The R dependence of the interactions shown in Figs, 
[l] and [2] is unlike the dipole-dipole 1/i? 3 interaction 
assumed in most theoretical treatments examining the 
many-body physics effects of these interactions. It is pos- 
sible to recover an interaction of this form by working at 
interatomic distances that are large compared to the crit- 



ical distance set by the condition where 25 = U(R C ), but 
it will be challenging to work in this regime with reason- 
able parameters. There are a number of requirements: 
the size of the interaction energy should be sufficiently 
large, in the 1-10 kHz range (similar to a BEC chemi- 
cal potential or the atom-atom interaction in a typical 
lattice); the spontaneous emission lifetime should be suf- 
ficiently long such that the many-body system can come 
into equilibrium (> 10 ms); the relevant interatomic dis- 
tance should match typical distances, either a lattice 
spacing or the average interparticle distance in a gas (0.5 
- 1 /im). It is difficult to achieve all these requirements for 
reasonable values of the coupling Rabi frequency. Ref [1 
suggests working at rather low principle quantum num- 
bers (n=20), but requires a large Rabi frequency (100 
MHz) and yields a short lifetime (5 ms). 

Producing standard dipolar physics with Rydberg- 
dressed atoms seems challenging at best. Another ap- 
proach is to look at the independent interactions that 
are possible with accessible experimental parameters and 
consider new many-body physics that is not possible for 
normal condensed matter systems. While such many- 
body states are beyond the scope of this work, we can 
speculate. If we consider the simpler case without any 
level crossings (Fig. [3|, we have an interaction that is re- 
pulsive and R- independent up to a critical radius, where 
it falls off rapidly towards zero. We can consider the im- 
pact of this interaction on a gas and on atoms trapped 
in an optical lattice. If the density was low enough for a 
gas of atoms, this interaction would look essentially like 
a contact interaction. At very high densities where the 
interatomic spacing is small compared to R c the inter- 
actions just provide an overall energy shift, presumably 
with little impact on the system. In an intermediate den- 
sity region, however, the cutoff at R c can be expected 
to have significant impact on the correlation functions 
of the gas, perhaps leading to exotic phases and order 
in the cloud. For an optical lattice, the addition of a 
nearest neighbor interaction to the Bose-Hubbard model 
leads to such phenomena as a supersolid phase [5 . With 
the Rydberg-dressed atom interactions we could have a 
situation where the first n-ncighbor interactions are the 
same, followed by little or no interactions. As this has 
no condensed matter analog, prediction is difficult, but 
novel phases seem likely to exist. Finally, because param- 
eters exist that achieve reasonable interaction strengths 
and lifetimes with dctunings < 1 GHz, we can also con- 
sider introducing spin dependence into such interactions. 
If a two-level spin system is established between two sub- 
levels in the different hyperfine states of an alkali atom 
such as Rb, we can map all these interactions onto a spin- 
dependent model, as only one of the levels would interact 
with the dressing laser. 

One important and simplifying assumption we have 
made with this treatment is ignoring any short range 
physics, i.e. assuming a single dipolar coupling at small 
R. This assumption will not be valid when the dipole- 
dipole shift is on the order of the spacing between neigh- 
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boring levels, and the energy spectrum becomes ex- 
tremely crowded and complicated at short range (see [18] 
for example). More detailed calculations beyond simple 
dispersion expansions of the dipole coupling are needed to 
understand the energy spectrum as well as potential loss 
mechanisms. The effect of the dressing at short distance 
may be more robust than it seems, however, because it 
depends solely on the blockade effect. As long as the 
interactions dominate the detuning their actual value is 
irrelevant. Of course if the coupling laser couples directly 
to a short-range molecular state, potentially large losses 
would ensue. 
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FIG. 5: The energy levels of the three-atom system in 
the one-dimensional configuration for £1/2tt — 10 MHz, 
S/2n = 100 MHz, and c 3 /2tt = -1000 MHzx^m 3 . 



THREE INTERACTING DRESSED ATOMS 



So far, we have examined only the simplest interact- 
ing system: two interacting dressed atoms. In order 
to extend these calculations to many-body systems in 
the future, it is necessary to explore a more compli- 
cated system of atoms. The next level of complextiy 
would be a system of three interacting dressed atoms. 
We can examine the eigenvalues for a system of three 
dressed atoms, where the Hamiltonian has multiple in- 
dependent interaction terms on the diagonal, in two sim- 
ple configurations: a one-dimensional configuration that 
examines next-nearest neighbor interactions, and a two- 
dimensional configuration. The one-dimensional config- 
uration is a line of the three atoms, with distance R 
between adjacent atoms, and the two-dimensional con- 
figuration is each atom on the vertex of an equilateral 
triangle with side R. 

The Hamiltonian of atoms arranged in the two- 
dimensional configuration has two independent diagonal 
terms. Because the atoms are equidistant, the diagonal 
terms for each amount of excitation are degenerate and 
we can reduce the 8x8 matrix to a 4 x 4 matrix 
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FIG. 4: The energy levels of the three-atom system in 
the two-dimensional configuration for U/2tt = 10 MHz, 
S/2n = 100 MHz, and c 3 /2tt = -1000 MHzx^m 3 . 
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using the basis: \ggg), 1/V3(\rgg) + \grg) + 
\ggr) , l/\/3(\rrg) + \grr) + \rgr),and\rrr). In this 
case, there are avoided crossings for 28 — —C3/R 3 , as 
well as for 36 = —Scs/R 3 for negative detunings and 
attractive interactions. For a given detuning, the second 
of these crossings occurs at a larger interparticle spacing, 
and is a narrower avoided crossing. This means that 
the larger interaction occurs at the first, closer, radius. 
Figure [4] shows the energy levels of this system. For 
negative C3, the interaction energy looks similar to the 
two-atom case, with an enhanced value for the same 
Rabi frequency and detuning. 

In the one-dimensional case, where the three atoms lie 
on a line, we find an additional avoided crossing, where 
2S = — C3/(8i? 3 ). This crossing is at an even closer ra- 
dius, and is wider than the other two. Figure [5] shows 
these avoided crossings as part of the energy level plot. 
In this reduced-dimension case, we see that the addition 
of a next nearest neighbor breaks the degeneracy among 
some of the diagonal terms, requiring the full 8x8 ma- 
trix for evaluation. Further expanding these calculations 
to consider a system of four atoms leads only to addi- 
tional avoided crossings at longer ranges in the energy 
spectrum. 



VI. TIME-DEPENDENT EFFECTS 

Although up until this point we have been concerned 
with the steady-state dressed atoms, since we have ex- 
pressions for the eigenvalues, it is simple to calculate 
the time dependence for a sudden turn-on of the laser 
fields. This can be done by projecting the ground state 
of the bare atom onto the dressed atom basis, evolving 
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the dressed-atom states in time with the phase factor 
exp(— itit) and then projecting back to the bare-atom 
basis. Figure [6] shows three cases: the evolution of the 
states \gg), l/y/2(\gr) + \rg)), and \rr) for a.) no atom- 
atom interactions, b.) atom-atom interactions with 6 = 
(single-atom resonance), and c.) atom-atom interac- 
tions on the two-atom resonance (26 = U(R)). For case 
a.) we see Rabi flopping of two independent atoms, with 
all the population in \rr) for a 7r-pulse, and population 
in all three states during the rest of the evolution. For 
case b.) the blockade effect is apparent, in that there is 
no observable excitation of the doubly excited state, and 
the system Rabi flops between \gg) and \gr) + \rg). For 
case c.) there is no population of single Rydberg atoms 
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FIG. 7: The population of the doubly-excited state after 
a pulse of duration t — 16/J7 at different detunings. 
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FIG. 6: The time-dependent evolution of the three 
eigenstates for atoms separated by 1 /jui for a.) no 
atom-atom interactions, b.) non-zero interations, but no 
detuning, and c.) both detuning and interactions. The 
colors show the evolution of the states \gg) (dashed), 
l/V2(\gr) + |rg))(solid), and |rr)(dotted). 



and the system flops between \gg) and \rr). Note the 
time scale is longer, because this is a two-photon cou- 
pling with a correspondingly smaller Rabi frequency. 

As can be seen in the figure, using a 7r/2-pulse with 
respect to the two-photon Rabi frequency should be a 
simple way to create an entangled state, \gg) + \rr) . Only 
atom pairs that are at the correct distance to be two- 
photon resonant will Rabi flop, so this should allow the 
creation of entanglement of pairs of atoms with ground- 
state atoms in between, for example. Figure [7] shows the 
population at the end of a pulse of duration t = 16/f2 
as a function of the detuning. For small detunings the 
blockade is evident as there is practically no population 
of the state \rr). At the two-photon resonance condition, 
all of the population can get excited into \rr) for the 
appropriate pulse duration. 



VII. CONCLUSION 

In summary we have considered creating Rydberg- 
dressed atoms, with small admixtures of Rydberg char- 
acter, to study dipolar many-body physics. We find that 
the blockade phenomena prevents large interaction ener- 
gies at small distances, set by a critical radius where the 
detuning is equal to the dipolar interaction energy. The 
resulting i?-dependent interaction energy can have suffi- 
cient size to be relevant for many-body physics of cold 
atoms, but its functional form is fundamentally differ- 
ent from the 1/i? 3 dipolar interaction that one expects. 
Instead it tends to a constant value at small R with a 
crossover to 1/R 3 at the critical distance. This form is 
unlike any interactions in condensed matter physics and 
may open up the possibility of novel many-body states. 
While these calculations point to experimental parame- 
ters that are challenging to achieve, in the future we hope 
to implement an experiment based on these calculations. 
A recent paper, [19], discusses further many-body effects 
of the interactions of dressed Rydberg atoms. 

We acknowledge discussions with J. Porto, S. Das 
Sarma, and G. Pupillo. This work is partially supported 
by the JQI NSF-PFC. 
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